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ABSTRACT 



Aims. To investigate the Ba and Eu abundances for a sample of 63 giant stars in the globular cluster M15. This is the largest sample of 
M15 giants stars for which Ba abundances have been determined and, due to the target selection of the original research programme, 
the Ba abundances are complete along the red giant branch. 

Methods. Stellar parameters were taken from the previous key study and a microturbulence-surface gravity relation was determined 
for precise measurement of the Ba line at 6496.898 A, which has a high sensitivity to microturbulence. Element abundances for Ba, 
La, Eu, Ca, Ni and Fe were calculated using spectrum synthesis and equivalent widths techniques. 

Results. A bimodal distribution in Ba, Eu and La abundances was found within the sample. The low Ba,Eu,La mode had mean 
abundances of <[Ba/H]>=-2.41+0.16, <[Eu/H]>=-1.80±0.08 and <[La/H]>=-2.19±0.13 while the high Ba,Eu,La mode had mean 
abundances of <[Ba/H]>=~2.00+0.16, <[Eu/H]>=-1.65±0.13 and <[La/H]>=-1.95±0.1 1. 

Conclusions. Both modes are indicative of a pollution scenario dominated by the r-process, hence contributions from explosive 
nucleosynthesis of massive stars. There may be evidence of further enhancement by another heavy element process and of potential 
anticorrelations in Na-O for both modes indicating a complex formation and evolution history for Ml 5. 
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1. Introduction 

The globular clusters of the Milky Way provide a wealth of 
information regarding the evolution of stars and the galaxy. 
Recent advances in technology and analysis techniques have re- 
vealed that many globular clusters show evidence of multiple 
stellar populations through photometry and chemical abundance 
distributions. A recent review of these discoveries is given in 
Gratton et al. (2012). Disentangling the different stellar popula- 
tions by their chemical signatures allows astronomers to identify 
the pollution and evolution mechanisms at work in these com- 
plex stellar systems. 

The origins of many of the light (Z < 30) and heavy (Z > 30) 
elements can be attributed to key stages in stellar evolution, 
whereby each stage produces unique chemical signatures, de- 
pending on the mass and metallicity of the progenitor star. The 
contributions from the different stages of stellar evolution are 
determined by identifying the quantity and distribution of these 
elements in the observed stellar population and considering then 
in combination with the knowledge of the sites of nucleosynthe- 
sis in which the elements are produced. 

The metal- poor globular cluster Ml 5, with [Fe/H]=- 
2.31+0.06 dex (icarretta et"al1 l2009b). has proven to be an in- 
teresting case for study due in particular to star-to-star variations 
in [Eu/H] and [Ba/H] that have been found. Eu is predominantly 



an r-process element, formed due to the rapid capture of neu- 
trons by a seed nuclei (Fe) in a high density neutron environ- 
ment. The key stellar evolution nucleosynthesis site for the r- 
process is explosive nu cleosynthesis in ma ssive stars, such as 
in Type II Supernovse (ISneden et al.l l2008l) . which is also the 
key production site for the a elements (O, Mg, Si, Ca, Ti ...) 
(Edvardsson et al. 1993). Hence comparing Eu and a elements 
may be used to probe the mass-range of the supemovae progeni- 
tors. Ba is a neutron capture element that can be produced by the 
r-process, but is predominantly produced in the solar system ma- 
terial by the i-process, which is the slow capture of neutrons by 
a seed nuclei. The typical nucleosynthesis site for the .s-process 
is within thermally pulsing asymptotic giant branch (AGB) stars 
jBusso etalj|2001h and it is possibly also significant in massive 
fast-rotating stars (flliadis 2007). Hence the comparison of s- and 
r-process element abundances also provides a means by which to 
discriminate between the contributions from these very different 
nucleosynthesis sites. 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile. Proposal 
ID: 080.B-0784(A) 



The study of lSneden et al.l (1 19971 hereafter, S97) investigated 
the abundances of both light (O, Na, Mg, Al, ...) and heavy 
(Ba, Eu) elements in a sample of 18 bright giants in M15. Anti- 
correlations typical to globular cluster were observed between 
Na and O, and between Al and Mg. Fe, Ca and Si were found to 
have uniform abundances for the sample. Eu and Ba were found 
to be strongly correlated with a possible bimodal distribution of 
both abundances, and to each have a large spread in values on 
the order of 0.6 dex. But no obvious correlation was found be- 
tween these heavy elements and the light elements. Hence it was 
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concluded that Ba and Eu had a primordial origin due to Type 
II Supernova of a stellar mass high enough to produce r-process 
elements but not sufficiently high to produce a elements, hence 
the lack of any correlation between those two types of elements. 

This work was then followed by two further studies, 
ISneden et all (l2000aL hereafter, SOOa) and Sneden et al. (2000b, 
hereafter, SOOb). In the former st udy three of the star s at the tip 
of the giant branch analysed in ISneden et al.l (11997 ) were re- 
observed to investigate further the neutron-capture abundances 
in these stars. The stars were chosen as having similar stel- 
lar parameters but v ery different Eu and Ba abundances as per 
ISneden et al.l (Il997h . These differences were confirmed and the 
heavy element signatures that were derived closely followed that 
produced primarily by the r-process. 

In the latter study a further 31 giants stars in Ml 5 were anal- 
ysed to investigate the proton-capture (Na), a-capture (Ca, Si, 
Ti) and neutron-capture (Ba) element abundances. It was found 
that there was a large range in Na abundances, which was at- 
tributed to proton-capture on Ne occurring internally to the ob- 
served st ars, as well a s a lar ge range in Ba comparable to that 
found in ISneden et al.l (11997 ) confirming that result and the bi- 
modal distribution of Ba. Eu abundances could not be measured 
for this sample. 

lOtsuki et al.l (l2006i hereafter, O06) investigated the heavy el- 
ement abundances of seven giant stars in M15 that had also been 
analysed in S97. They found Fe and Eu abundances in agree- 
ment with S97 and extended the analysis to Zr, Y and La. They 
found that La correlated with Eu, as expected, but observed an 
anti-correlation of Zr and Y with Eu. From this they concluded 
that while the La (and Ba) correlation with Eu implies a pri- 
mary contiibution of neutron-capture elements by an r-process 
only mechanism (e.g. explosive nucleosynthesis within massive 
stars) this cannot also explain the anticorrelation of Zr and Y 
with Eu. They argue that some other process enriched the inter- 
stellar medium in Ml 5 in light neutron-capture elements before 
the enrichment of the heavy neutron-capture elements by super- 
nova. 

iPreston et al.l (l2006i hereafter, P06) made an investigation of 
six red horizontal branch (RHB) stars in M15. Both light and 
heavy element abundances were determined, for the M15 stars 
similar results regarding the evidence of star-to-star scatter in the 
heavy element abundances were found. Also noted was that the 
RHB stars of M15 were found to be ~0.2 dex more metal-poor 
than red giant branch (RGB) stars in the same cluster 

Most recently, a detailed investigation into neutron-capture 
element abundances in giant stars in M15 was carried out by 
ISobeck et all (1201 1. hereafter, Skll). The three RGB stars of 
SOOa and the six RHB stars of P06 were re-analysed in a con- 
sistent manner in order to investigate the nature of the heavy 
element abundances in M15 giants stars, and to investigate the 
metallicity offset between the RHB and RGB as observed in 
P06. The abundances for 40 species of both light and heavy el- 
ements were deiived for these nine stars. As had been observed 
previously the abundances of proton-capture elements showed 
no correlation with the abundances of the neutron capture ele- 
ments. In particular a metallicity offset between the RHB and 
RGB was found, and the large spread in Ba and Eu was con- 
firmed. Ho wever the bimodal di stribution of Ba and Eu that was 
proposed in lSneden et al.l (Il997h was not observed. The analysis 
concluded that the heavy neutron-capture element abundances 
agreed well with a scaled-solar r-process distiibution, but the 
light neutron-capture element abundances did not agree with ei- 
ther a scaled-solar s- or r-process distribution. This was therefore 



evidence of some other nucleosynthetic process, under the um- 
brella of the Lighter Element Primary Process (LEPP), at work. 

The largest sample of M15 stars ana lysed spectroscopically 
to-date was carried out bv lCarretta et al.l (i2009b, hereafter, C09) 
as part of an analysis of light element abundances for stel- 
lar samples within 15 galactic globular clusters observed using 
FLAMES/GIRAFFE. Fe, Na and O abundances were investi- 
gated for 84 stars in M15. In total the research programme ob- 
tained homogeneously determined abundances for Na and O for 
1235 stars in 19 globular clusters. An anti-correlation between 
Na and O was observed for all the clusters, including M15. This 
provided evidence of a pristine first generation of stars (~30% 
of the observed cluster stars) and a second generation that had 
formed from material enriched in H-burning pro ducts by now- 
extinct massive stars of the first stellar generation. ICarretta et al.l 
(2009a) followed up this study with high resolution observations 
of 200 red giants across 17 GCs, for which 13 stars were ob- 
served in M15 from the C09 sample. 

Of particular r elevance is the study undertaken in 
iD'Orazi etaP (1201 0'. hereafter, DOlO) which investigated Ba 
abundances in GC stars directly from the observations of C09. 
This study measured Ba for 1200 stars across the 15 GCs, in- 
cluding measurements for 57 stars in the M15 sample. M15 was 
the only sample for which they found a large variation in the 
Ba abundances between the stars confirming the large variation 
found in S97 and SOOb. 

The study carried out here seeks to extend the sample of gi- 
ant stars in M15 for which heavy neutron-capture element abun- 
dances have been determined and consider their distribution in 
M15 in light of the previous studies. 

Section |2] describes the observations and presents the stellar 
parameters for the dataset under analysis. Section [3] describes 
the determination of the stellar parameters as per Carretta et alj 
(2009b) and the treatment of the microtubulence, ^, in this study. 
Section |4] describes the determination of the chemical abun- 
dances for this sample. Section|5]present the results of the chem- 
ical abundance analysis. Section |6] discusses possible scenarios 
based on the results. Section |7]concludes the paper. 



2. Observations 

Observations of 63 red giant branch (RGB) stars in M15 were 
taken using the HR14A setup of GIRAFFE on the VLT in 
November 2009. The spectra were wavelength calibrated as 
from the ESO pipeline and sky-subtraction was performed 
using IRAF. Normalisation and radial velocity determination 
were performed using tools developed for the AMBRE Project 
(Worley et al. 2012). The signal-to-noise (S/N) ranged from 50 
to 170 across the sample. 

These targets were taken from C09 in which the stellar pa- 
rameters of effective temperature (7"gff), surface gravity (logg) 
and metallicity ([Fe/H]) as well as individual chemical abun- 
dances for Na and O were obtained for a total sample of 84 stars 
in M15. Table[T]lists the star ID and stellar parameters of the 63 
stars from C09 that were analysed here as taken from Table 3 in 
C09. 

The DOlO sample was also a subsample of C09. Of the 57 
DOlO stars, 38 were also observed for this study. DOlO could 
not analyse the complete C09 sample for Ba abundances as not 
all the stars were observed using the HR13 grating, in which the 
Ba spectral feature at 6141.73 A is present. 

The colour-magnitude diagram (CMD) and Hertsprung- 
Russell (HR) diagram of the 63 stars from C09 that have been 
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Table 1. The star ID, effective temperature, surface gravity, model metallicity, microturbulence, radial velocity, abundances derived 
from Fe I and Fe II lines, Na and O abundances and their associated root mean square error (rms) for each of the Ml 5 stars analysed 
here as determined in C09 



ID 


^eff 
(K) 


log^f 
(dex) 


[A/H] 
(dex) 


(kms-') 


(kms-') 


[Fe I/H] 
(dex) 


rms 


[Fe II/H] 
(dex) 


rms 


[0/Fe] 
(dex) 


rms 


[Na/Fe] 
(dex) 


rms 


40825 


4313 


0.65 


-2.33 


2.25 


-113.91 


-2.33 


0.14 


-2.34 


- 


0.24 


0.11 


0.44 


0.07 


4099 


4324 


0.69 


-2.32 


2.03 


-109.43 


-2.32 


0.11 






0.54* 


0.13 


-0.07* 


0.01 


43788 


4325 


0.52 


-2.36 


2.40 


-96.93 


-2.36 


0.10 


-2.46 


0.08 


0.39 


0.02 


0.13 


0.06 


31914 


4338 


0.58 


-2.20 


2.09 


-112.77 


-2.20 


0.10 


-2.21 


- 


-0.32 


- 






41287 


4373 


0.78 


-2.36 


1.10 


-116.25 


-2.36 


0.07 










0.39 


0.01 


34519 


4416 


0.87 


-2.38 


1.62 


-98.04 


-2.38 


0.11 


-2.46 


- 


0.46 


0.04 


0.11 


0.08 


37215 


4445 


0.98 


-2.31 


1.28 


-111.36 


-2.31 


0.01 










0.73 


- 


34995 


4468 


1.00 


-2.34 


0.79 


-115.14 


-2.34 


0.08 










0.64 


- 


3137 


4486 


1.08 


-2.35 


2.27 


-109.50 


-2.35 


0.13 


-2.29 


- 


0.29 


0.08 


0.58 


0.04 


42262 


4528 


1.13 


-2.33 


2.50 


-104.50 


-2.33 


0.08 










0.90 


0.13 


26751 


4533 


1.18 


-2.44 


1.04 


-98.72 


-2.44 


0.08 










0.55 


0.06 


38678 


4554 


1.23 


-2.35 


1.28 


-109.01 


-2.35 


0.08 


-2.46 


0.01 


0.13 


- 


0.56 


0.03 


2792 


4567 


1.26 


-2.32 


1.50 


-100.93 


-2.32 


0.13 






0.32* 


0.00 


0.20* 


0.06 


39752 


4598 


1.31 


-2.36 


1.30 


-108.38 


-2.36 


0.13 


-2.39 


0.03 


0.26 


- 


0.29 


0.15 


34456 


4621 


1.37 


-2.34 


1.40 


-97.02 


-2.34 


0.20 






0.22 


- 


0.34 


0.01 


34961 


4627 


1.41 


-2.32 


1.85 


-113.28 


-2.32 


0.15 


-2.36 


0.02 


-0.01 


- 


0.65 


0.03 


39787 


4630 


1.35 


-2.43 


2.10 


-111.94 


-2.43 


0.13 






0.48* 


0.00 


0.13* 


0.03 


38329 


4664 


1.44 


-2.37 


0.86 


-109.04 


-2.37 


0.02 










0.16 


- 


45062 


4700 


1.54 


-2.31 


1.66 


-110.27 


-2.31 


0.11 


-2.23 


- 


0.03 


- 


0.63 


0.05 


20498 


4717 


1.61 


-2.49 


1.12 


-106.58 


-2.49 


0.08 


-2.44 


- 


0.25 


- 


0.54 


0.04 


21948 


4746 


1.66 


-2.49 


0.20 


-108.68 


-2.49 


0.05 










0.46 


- 


28510 


4754 


1.64 


-2.42 


0.14 


-100.52 


-2.42 


0.02 






0.23* 


0.00 


0.07* 


0.07 


29264 


4810 


1.77 


-2.35 


0.84 


-101.86 


-2.35 


0.14 






0.40* 


0.00 


0.01* 


0.00 


18815 


4832 


1.83 


-2.29 


0.80 


-108.68 


-2.29 


0.10 






0.66* 


0.00 


0.00* 


0.00 


42362 


4952 


2.03 


-2.41 


1.00 


-110.67 


-2.41 


0.09 






0.27 


- 


0.38 


0.17 


29480 


4968 


2.10 


-2.33 


1.16 


-104.60 


-2.33 


0.16 






0.30 


- 


0.45 


0.14 


37931 


4976 


2.10 


-2.37 


1.21 


-111.25 


-2.37 


0.05 










0.27 


0.08 


31791 


4978 


2.13 


-2.28 


0.10 


-104.92 


-2.28 


0.03 










0.14 


- 


33484 


5036 


2.22 


-2.44 


1.20 


-109.75 


-2.44 


0.17 










0.14 


- 


10329 


5060 


2.25 


-2.38 


1.25 


-109.91 


-2.38 


0.11 






0.42 


- 


0.53 


- 


27889 


5063 


2.28 


-2.32 


1.15 


-101.14 


-2.32 


0.13 






0.63 


- 


0.04 


0.11 


39741 


5076 


2.31 


-2.25 


1.42 


-112.10 


-2.25 


0.20 






0.43 


- 


0.38 


0.18 


42674 


5076 


2.29 


-2.43 


1.35 


-96.28 


-2.43 


0.06 










0.67 


- 


17458 


5081 


2.33 


-2.31 


0.92 


-110.56 


-2.31 


0.19 






0.44 


- 


0.16 


0.05 


40086 


5087 


2.36 


-2.32 


1.48 


-122.86 


-2.32 


0.10 






0.41 


- 


0.37 


0.07 


36274 


5105 


2.35 


-2.35 


1.16 


-100.30 


-2.35 


0.13 






0.50 


- 


0.34 


0.04 


26759 


5106 


2.34 


-2.35 


1.16 


-101.07 


-2.35 


0.20 










0.24 


- 


29521 


5129 


2.46 


-2.30 


0.03 


-101.45 


-2.30 


0.02 










0.19 


0.08 


41987 


5130 


2.41 


-2.33 


2.20 


-108.07 


-2.33 


0.09 










0.39 


- 


1939 


5149 


2.48 


-2.41 


0.88 


-101.98 


-2.41 


0.03 






0.32 


- 






34272 


5151 


2.47 


-2.25 


1.24 


-108.35 


-2.25 


0.10 






0.47 


- 


0.38 


0.13 


35961 


5156 


2.44 


-2.41 


0.40 


-107.90 


-2.41 


0.10 










0.38 


0.06 


31010 


5164 


2.46 


-2.32 


0.66 


-115.45 


-2.32 


0.08 










0.61 


0.03 


23216 


5216 


2.57 


-2.32 


0.14 


-104.01 


-2.32 


0.02 














18240 


5229 


2.59 


-2.30 


0.94 


-103.56 


-2.30 


0.03 










0.56 


0.20 


32942 


5236 


2.61 


-2.25 


1.82 


-108.90 


-2.25 


0.10 










0.38 


0.04 


32979 


5236 


2.57 


-2.35 


1.97 


-106.56 


-2.35 


0.04 










0.17 


- 


2411 


5239 


2.65 


-2.29 


1.07 


-110.01 


-2.29 


0.22 














23153 


5242 


2.59 


-2.26 


1.61 


-108.68 


-2.26 


0.10 














8927 


5250 


2.64 


-2.42 


0.16 


-104.84 


-2.42 


0.08 










0.30 


0.10 


18770 


5260 


2.63 


-2.38 


1.95 


-103.56 


-2.38 


0.10 














22441 


5270 


2.70 


-2.35 


0.35 


-108.68 


-2.35 


0.09 














18508 


5276 


2.68 


-2.26 


2.61 


-103.56 


-2.26 


0.09 














9608 


5291 


2.68 


-2.33 


0.31 


-104.47 


-2.33 


0.12 






0.41 








29436 


5312 


2.81 


-2.41 


1.39 


-98.38 


-2.41 


0.11 














4989 


5319 


2.72 


-2.38 


0.10 


-104.70 


-2.38 


0.06 














3825 


5323 


2.76 


-2.26 


0.72 


-109.12 


-2.26 


0.12 


-2.24 


0.34 










40762 


5325 


2.78 


-2.34 


1.17 


-115.00 


-2.34 


0.07 


-2.21 








0.66 




8329 


5344 


2.86 


-2.41 


1.77 


-106.98 


-2.41 


0.02 














2463 


5349 


2.84 


-2.31 


1.01 


-105.89 


-2.31 


0.12 


-2.27 




0.48 




0.40 


0.08 


19346 


5349 


2.80 


-2.34 


2.70 


-105.99 


-2.34 


0.04 










0.94 


0.03 


28026 


5362 


2.81 


-2.38 


1.80 


-111.48 


-2.38 


0.15 














31796 


5374 


2.87 


-2.25 


1.80 


-110.52 


-2.25 


0.08 










0.08 





* [O/Fe] Sl [Na/Fe] taken from Carretta et al. (2009a' 
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Fig. 1. The CMD and HR diagram of the key samples of giant stars in M15 as per the legend. The colours and parameters for 
W13 and DOlO both come from C09. 



observed and analysed here are shown in Figure [T] The locations 
of stars from previous key studies of M15 are also shown for 
comparison including the DOlO subsample of C09. 

The wavelength range of the HR14 grating covers from 
~6308 A to ~6700 A, which includes the key neutron-capture 
features of Ba II at 6496.898 A, La II at 6390.4770 A, and Eu II 
at 6437.640 A and 6645.127 A. At this resolution of R~14,500 
these three features are observable although the Eu II features are 
weak. The astrophysical quantities associated with these lines, as 
well as those for the other lines measured in this study, are listed 
in Table |2] Table [3] lists the hyperfine structure decomposition 
for the Ba II, La II and Eu II spectral lines as well as the ref- 
erences from which the atomic information was obtained upon 
which the calculations were based. 



3. Stellar Parameters 

The stellar parameters for each of the stars observed here had 
been determined in C09. The T^ff was determined in a two step 
process whereby a relation between T^f^ and the V (or K) mag- 
nitude was developed for th e cluster using V - K photometry 
and the Alonso et al ] (119991) calibration. This relation was then 
used to calculate the T^f^ for each star in the sample. The logg 
was determined using Tgff , the apparent magni tude, the distanc e 
modulus and bolometric corrections from Alonso„et al.J (Il999h . 
The ^ was determined by eliminating any trend of the abundance 
of iron with expected line strength based on the equivalent width 
measurement of Fe I lines. See C09 and references therein for 
more detailed descriptions of the methods employed. 

However, as can be noted from Table [1] the resulting micro- 
turbulence velocities for the C09 M15 sample range from 0.03 to 
2.7 kms"' which is a very wide and rather implausible range for 
this sample of giant stars. This wide range most probably results 
from too few and too weak Fe lines being available in the spec- 
tral range used in the C09 study for this very metal-poor cluster 
These ^ values anyway have almost no impact on the resulting 
abundances of Fe, O and Na investigated in C09, since the lines 
analysed were on the linear part of the curve of growth. On the 
other hand, the Ba line under analysis here is strong enough for 
the microturbulence to play a role in desaturating the lines. 



Table 2. List of Fe, Ca, Ni, Ba, Eu and La lines used in this 
analysis with the associated excitation potential (x), log gf and 
method of analysis used either spectrum synthesis (SS) or equiv- 
alent width curve of growth (EW). The reference solar abun- 
dances are also listed. 



Wavelength (A) 


Species 






Reference 


Method 


Sol* 


6439.0700 


20.0 


2.526 


0.470 


NIST 


SS 


6.41 


6471.6600 


20.0 


2.526 


-0.590 


NIST 


SS 




6318.0171 


26.0 


2.453 


-1.804 


NIST 


EW 


7.54 


6322.6851 


26.0 


2.588 


-2.426 


NIST 


EW 




6335.3304 


26.0 


2.198 


-2.177 


NIST 


EW 




6336.8239 


26.0 


3.686 


-0.856 


NIST 


EW 




6344.1487 


26.0 


2.433 


-2.923 


NIST 


EW 




6355.0286 


26.0 


2.845 


-2.291 


NIST 


EW 




6393.6009 


26.0 


2.433 


-1.576 


NIST 


EW 




6408.0179 


26.0 


3.686 


-1.018 


NIST 


EW 




6411.6489 


26.0 


3.654 


-0.718 


NIST 


EW 




6421.3504 


26.0 


2.279 


-2.027 


NIST 


EW 




6430.8460 


26.0 


2.176 


-2.006 


NIST 


EW 




6475.6240 


26.0 


2.559 


-2.942 


NIST 


EW 




6481.8699 


26.0 


2.279 


-2.984 


NIST 


EW 




6518.3667 


26.0 


2.832 


-2.298 


NIST 


EW 




6575.0154 


26.0 


2.588 


-2.710 


NIST 


EW 




6592.9134 


26.0 


2.728 


-1.473 


NIST 


EW 




6593.8701 


26.0 


2.433 


-2.422 


NIST 


EW 




6432.6800 


26.1 


2.891 


-3.500 


NIST 


EW 




6456.3760 


26.1 


3.903 


-2.190 


NIST 


EW 




6516.0530 


26.1 


2.891 


-3.372 


NIST 


EW 




6643.6300 


28.0 


1.675 


-2.300 


NIST 


SS 


6.29 


6496.8980 


Ball 


0.604 


-0.407 


012° 


SS 


2.25 


6390.4770 


La II 


0.321 


-1.410 


LOla* 


SS 


1.25 


6437.6400 


EuII 


1.320 


-0.320 


LOlb" 


SS 


0.60 


6645.1270 


EuII 


1.380 


0.120 


LOlb" 


SS 





* Solar reference abundances taken from Lodders ( 20Q3I) 

a) IDavid50n etai]fT993) : bl lLawleretaLgOOla) : c) ILawler et aHDOOlbf) 



The sensitivity of the 6496.898 A Ba line to ^ is clearly 
shown in the error analysis in Table |5] In a preliminary analysis 
using the ^ values determined in C09, a clear trend of Ba abun- 
dance with ^ was found. This is a typical problem when dealing 
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Table 3. Hyperfine structure decom positi on based on atomic informatio n obtained from the associated references for the Ba II 
(lRuttenlll978l:lDavidson et al.lll992h . La II ( iLawleret alJl2001ah and Eu II (lLawleretal.ll2001bl) spectral lines used in the spectiTim 
synthesis of the M15 giant stars. 



Wavelength(A) 



Strengtii 



Wavelength(A) 



Strengtii 



La II ^6496 

6390.455.... 
6390.468.... 
6390.468.... 
6390.479.... 
6390.479.... 
6390.480.... 
6390.489.... 
6390.489.... 
6390.490.... 
6390.496.... 
6390.497.... 
6390.498.... 
6390.502.... 
6390.503.... 
6390.506.... 



;ir=0.321 eV, logg/=-1.410 



Ba II /16496.898: x=0-604 eV, log g/= 



-0.407 



-2.012 
-2.183 
-2.752 
-2.570 
-3.752 
-2.390 
-2.536 
-3.334 
-2.661 
-3.100 
-2.595 
-3.079 
-2.954 
-2.778 
-2.857 



'Ba: 
6496.898.. 
"5Ba: 
6496.887.. 
6496.890.. 
6496.895.. 
6496.903.. 
6496.905.. 
6496.908.. 
"SBa: 
6496.898.. 
'"Ba: 
6496.886.. 
6496.889.. 
6496.894.. 
6496.904.. 
6496.906.. 
6496.909.. 
"^Ba: 
6496.898.. 



-0.406 

-1.911 
-1.212 
-0.765 
-1.610 
-1.212 
-1.212 

-0.406 

-1.911 
-1.212 
-0.765 
-1.610 
-1.212 
-1.212 

-0.406 



Eu II ^6437.640: ;r= 1-320 eV, logg/=-0.320 



^Eu: 
6437.611.. 
6437.619.. 
6437.629.. 
6437.636.. 
6437.643.. 
6437.650.. 
6437.656.. 
6437.662.. 
6437.666.. 
6437.672.. 
6437.676.. 
6437.679.. 
6437.684.. 
6437.687.. 
6437.689.. 
6437.692.. 
'"Eu: 
6437.610.. 
6437.613.. 
6437.624.. 
6437.627.. 
6437.630.. 
6437.634.. 
6437.637.. 
6437.639.. 
6437.640.. 
6437.641.. 
6437.642.. 
6437.642.. 
6437.643.. 
6437.644.. 
6437.644.. 
6437.645.. 



Eu II .16645.127: ;i^= 1.380 eV, logg/=+0.120 



Eu: 



-0.960 


6645.071 


-0.517 


-2.191 


6645.078 


-1.823 


-2.191 


6645.085 


-3.467 


-1.070 


6645.097 


-0.593 


-1.998 


6645.104 


-1.628 


-1.998 


6645.112 


-3.149 


-1.181 


6645.118 


-0.672 


-1.956 


6645.126 


-1.583 


-1.956 


6645.133 


-3.077 


-1.287 


6645.137 


-0.754 


-2.010 


6645.144 


-1.635 


-2.010 


6645.150 


-3.245 


-1.377 


6645.152 


-0.839 


-2.207 


6645.158 


-1.830 


-2.207 


6645.163 


-0.921 


-1.428 


'"Eu: 






6645.072 


-1.823 


-0.960 


6645.074 


-0.517 


-2.191 


6645.074 


-3.467 


-2.191 


6645.087 


-0.593 


-1.070 


6645.089 


-1.628 


-1.998 


6645.094 


-3.149 


-1.998 


6645.096 


-0.672 


-1.181 


6645.101 


-1.583 


-1.956 


6645.104 


-0.754 


-1.956 


6645.107 


-3.077 


-2.010 


6645.109 


-0.839 


-1.287 


6645.110 


-1.635 


-2.207 


6645.113 


-0.921 


-1.428 


6645.116 


-1.830 


-1.377 


6645.116 


-3.245 


-2.010 






-2.207 







with the measurement of strong lines such as those of Ba in gi- 
ants. In SOOb a constant ^ of 2.0 kms"' was assumed while Ski 1 
treated ^ so as to eliminate any trend for the strong lines. 

For this sample we therefore re-determined ^ using the avail- 
able Fe I lines as listed in Table |2] Hence the equivalent widths 



for Fe I and Fe II lines were measured for those stars in which a 
sufficient number of unblended lines could be identified. For the 
majority of stars the combination of low S/N and high tempera- 
tures meant a insufficient sample of Fe lines could be identified. 
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Table 4. Stellar Parameters from C09 for 1 3 stars for which the equivalent widths of Fe I and Fe II lines were measured and/or were 
analysed in previous studies as listed. The mean and standard deviation and the number of lines measured for Fe I and Fe II are 
listed. The standard deviation is given only for measurements of N>2. 



Star ID 


Other ID 


^eff 




[A/H] 






[Fe I/H] 


cr 


N* 


[Fe II/H] 


cr 


N 


Other Studies 


40825 


K386 


4313 


0.65 


-2.33 


2.25 


2.25 


-2.40 


0.10 


35 (22) 


-2.37 


_ 


2 


S97, SOOb, Ot06 


4099 


K341 


4324 


0.69 


-2.32 


2.03 


2.30 


-2.40 


0.11 


17 


-2.34 


_ 


2 


S97, SOOa, SOOb, Skll 


43788 


K825 


4325 


0.52 


-2.36 


2.40 


2.25 


-2.32 


0.12 


17 


-2.37 


_ 


2 


S97, SOOb 


31914 


K757 


4338 


0.58 


-2.20 


2.90 


2.10 


-2.27 


0.11 


17 


-2.33 




2 


S97, SOOb 


41287 


K702 


4373 


0.78 


-2.36 


1.10 


2.10 


-2.41 


0.13 


17 


-2.37 




2 


S97 


34519 




4416 


0.87 


-2.38 


1.62 


2.00 


-2.44 


0.10 


23 (9) 


-2.32 




2 




34995 




4468 


1.00 


-2.34 


0.79 


2.10 


-2.41 


0.13 


15 


-2.40 




2 




3137 


K387 


4486 


1.08 


-2.35 


2.27 


2.00 


-2.38 


0.12 


16 


-2.32 




2 


S97 


26751 


K709 


4533 


1.18 


-2.44 


1.40 


2.02 














SOOb 


38678 


K288 


4554 


1.23 


-2.35 


1.28 


2.00 


-2.47 


0.11 


20 (9) 


-2.38 




2 


SOOb 


2792 


K255 


4567 


1.26 


-2.32 


1.50 


2.20 


-2.47 


0.13 


11 


-2.37 


0.14 


3 


SOOb 


34961 




4627 


1.41 


-2.32 


1.85 


2.00 


-2.46 


0.13 


22 (14) 


-2.46 


0.09 


3 




38329 




4664 


1.44 


-2.37 


0.86 


1.80 


-2.41 


0.12 


11 


-2.28 




2 





# The numbers in brackets are the number of equivalent widths obtained from C09 



-2.1 
-2.2 
5 -2.3 
t-2.4 
-2.5 
-2.6 



-2.1 ^ 
-2.2 
-2.3 
-2.4 
-2.5 
-2.6 - 
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Fig. 2. Comparison of [Fe I/H] and [Fe II/H] values from C09 Fig. 3. Comparison of C09 ^ and logg values with the ^-\ogg 

and the values determined here by equivalent width measure- relation from K09, the ^ values determined here for the sample 

ment. The C09 [Fe I/H] are those used as input for the stellar of 12 stars and the resulting ^-log § relation that was derived, 
atmosphere models. 



Table |4] lists the C09 parameters for 13 stars of which 12 
had a sufficient sample of Fe lines that could be used to deter- 
mine the microturbulence. The microturbulence determined in 
this analysis {^wn) with the mean [Fe I/H], mean [Fe II/H], asso- 
ciated standard deviations and number of lines measured are also 
listed. For four of the spectra the equivalent widths of Fe I lines 
measured in C09 for the corresponding stars were also used and 
the number of lines for each are also indicated. Hence for each 
of these 12 stars the equivalent widths of all possible Fe I and 
Fe II lines were measured and abundances determined using the 
curve of growth analysis routine, abfind, in the one-dimensional, 
local thermodynamic equilibrium (LTE) radiative transfer code, 
MOOG (Sneden 1973). Stellar models for each star were gener- 
ated at the C09 stellar parameters using the Kurucz 1995 stellar 
model grid and the ATLASSQ model atmosphere calculation pro- 
gramme. 



' http ://kurucz. harvard. edu/programs .html 



Eight of this sample of stars had also been analysed in pre- 
vious studies, as well a ninth star for which Fe I lines could not 
be measured here. The papers associated with each star are also 
listed in Table g] 

Figure |2] compares the mean [Fe I/H] and mean [Fe II/H] 
values measured here with those measured in C09 for the 12 
stars, although in C09 Fe II lines could only be measured for 6 of 
the current sample. The errorbars on the W13 measurements of 
Fe II lines come from the uncertainty on the measurement of the 
Fe I line in the corresponding star, as the dispersion around the 
Fe I lines is a fair representation of the observational scatter for 
the measurement of weak lines in these spectra. Thus Figure |2] 
shows that there is agreement within 1 cr between C09 and W13 
for [Fe I/H] and [Fell/H]. 

Based on the ^wi'i derived for these 12 stars a f-log g relation 
was determined as the linear best fit to the 12 data points. This 
relation was then used to provide a ^wn value for the remaining 
stars in the sample calculated from their log g value. Figure |3] 
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shows the ^ and logg values as determined in C09 and in this 
study (W13) with the ^-\ogg relation and uncertainty limits de- 
termined in Kirby et al. (2009) (K09). The linear best fit to the 
^wi3 values is also shown, extended across the range of logg 
values for this sample. The equations for the K09 relation and 
relation derived here are as follows: 

= (2.13 + 0.05) - (0.23 + 0.03) logg, (1) 
= 2.386 -0.3067 logg. (2) 

This sample is tailored specifically to these mono-metallicity 
stars while the K09 relation was derived from a sample of stars 
covering a range of metallicities. However the relations are in 
reasonably good agreement and in the expected range for giant 
stars. 

4. Determination of Stellar Abundances 

As mentioned above the iron abundances based on Fe I and Fe II 
lines were determined using equivalent widths and curve-of- 
growth analysis in MOOG. For the determination of the abun- 
dances for Ca, Ni, Ba, La and Eu the techniques of spectrum 
synthesis were employed, again using MOOG and the Kurucz 
1995 stellar models interpolated to the C09 parameters but now 
with the ^wi3 values. The most recent laboratory log g/s were 
obtained, as listed in Table |2] and the hyperfine structure for the 
Ba II, La II and Eu II were taken into account as per the decom- 
position in Table[3]bas ed on the most recent atomic info rmation 
(lRutteniri 978; Davids on et al.lll992l:lLawler et al.ll2001allbh . 

4. 1 . Comparison Stars 

As shown in Table|4]there are nine stars in this sample that have 
been analysed in previous studies. Figure |4] shows the range in 
values obtained for these stars for their stellar parameters and 
individual element abundances. 

The stellar parameters of T^ff are mostly within 50 K in 
terms of agreement between the studies (Figure 2^). However 
for four of the stars the spread in log g is on the order of 
0.5 dex (FigurellJ?) and ^ is particularly spread out (~0.3 kms"') 
between the studies for the high luminosity, low logg stars 
(Figure Ht), which is reflective of the difliculty of deriving pa- 
rameters for these types of stars. These differences reflect the de- 
gree of variation in stellar parameters that can be derived for the 
same star using diff'erent observations, measurements and anal- 
ysis techniques. For example the three hottest stars have param- 
eters from C09 and SOOb only. For SOOb, if the stars were not 
from S97, then the parameters were inferred from photometry - 
T^ff, log g relations. Hence the C09 parameters and SOOb param- 
eters are in good agreement for Tggp and log g, compared to the 
next two coolest stars which have C09 and S97 parameters. S97 
derived parameters spectroscopically and while the T^f^ agree 
reasonably well with C09 the log g do not. However the good 
agreement in [Fe/H] between the studies is seen in Figure |4j;. 

Figure |4jl to h compare the abundances for Ca, Ni, Ba and 
Eu derived here with the abundances from the other studies for 
the comparison stars. There is reasonable agreement between 
the studies. Ca in particular shows good agreement between the 
studies at a reasonably constant absolute abundance. Ni is con- 
stant in abundance for this study but the three comparisons to 
two other studies show some variation. The range in abundance 
values between the stars for Ba and Eu is obvious already for 
this small sample of stars. For the chemical abundances the dis- 
agreement between studies is greatest for the high luminosity. 



reflecting the uncertainties in stellar parameters and chemical 
abundance determination for stars of low log g. There may also 
be uncertainties in the model atmospheres (eg. sphericity effects) 
and in non-local thermodynamic equilibrium effects associated 
with cool giants of low gravity. 

Table 5. Variations in measured abundances with changes in 
stellar atmospheric parameters for the sample of Fe I and Fe II 
lines, by equivalent width measurement, and the other key lines 
used in this study. 





AT'eff 
+50 K 


+0.5 dex 


-0.25 kms-' 


A[Fe I/H] 


0.10 


-0.06 


0.06 


A[Fe II/H] 


-0.02 


0.18 


0.01 


A[Ca I/Fe](6439.07) 


0.10 


-0.08 


0.14 


A[CaI/Fe](6471.66) 


0.06 


-0.06 


0.02 


A[Ni I/Fe] 


0.10 


-0.03 


0.07 


A[Ba II/Fe] 


0.05 


0.17 


0.23 


A[Eu II/Fe](6437.64) 


0.00 


0.16 


0.08 


A[Eu II/Fe](6645.13) 


0.00 


0.17 


0.00 


A[La II/Fe] 


0.03 


0.18 


0.01 



Table |5] gives the expected variations in the derived 
abundances for a Arg^=50 K, Alogg=+0.5 dex and A^-- 
0.25 kms"' for the light and heavy element lines being anal- 
ysed in this study. The neutral species show greater sensitivity 
to changes in Tgff, while the ionised species show greater sensi- 
tivity to changes in logg as expected. The Ca I line at 6471.66 A 
and the Ba II line both show a sensitivity to changes in ^, in par- 
ticular the Ba II line which varies by 0.23 dex, as is expected for 
strong lines. 

The possible effects of departures from LTE on the Ba 
spe ctral feature at 649 6 A h ave b een previously examine d 
by IShort & Haus childtl (l2006l) and lAndrievskv etaP (l2009l) . 
Short & Hauschildil2006|) analysed the non-LTE sensitivity of 
the Ba II 6496 transition for a variety of different stellar at- 
mospheric model and line transfer combinations in a typical 
halo giant (with T^ff - 4800 K and logg = 1.5). By and large, 
they found a metallicity dependence of the non-LTE effect with 
line strength (explicitly, a non-LTE weakening of transitions 
for the most metal-deficient models and conversely, a non-LTE 
strengthening of lines for the least metal-poor models). Yet, in 
particular for the 6496 transition, when the data for full non-LTE 
model/transfer are compared to those for LTE model/transfer, the 
magnitude of the change was relatively small on the resultant Ba 
abundance (with | logW^(NLTEFull)- logW,,(LTE)| < 0.03). 

On the other hand, Andrievskv et al. (2009) employed a sin- 
gle LTE atmospheric model with non-LTE radiative transfer 
to study the non-LTE effect on Ba II lines as a function of 
both metallicity and temperature in extre mely metal-poor gi- 
ants. Their results mostly confirm those of IShort & Hauschildtl 
(2006) (a correlation between metallicity and non-LTE depar- 
ture was established with the largest effect seen for warm, metal- 
poor stars). Specifically for the Ba II 6496 A feature, they found 
a maximum negative correction factor (NLTE-LTE) of ^-0.1 
for a giant with [Fe/H] =-2.0 and T^ff - 4700 K and a max- 
imum positive (NLTE-LTE) correctio n of ^-0.13 for a gian t 
with [Fe/H ]=-2.5 and Teff = 5500 K. I Andrievskv et al.1 (|2009|) 
note, however, that the correction factor cannot be used directly 
to modify the Ba abundance (and rather, the comparison of non- 
LTE to LTE line profiles is requisite). Consequently, in consid- 
eration of these findings as well as those of the current study. 
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Fig. 4. Stellar parameters and chemical abundances against V magnitude for the 9 comparison stars in the M15 sample. The values 
from C09, W13, S97, SOOa, SOOb, O06 and Ski 1 ai-e shown: ajT^ff, b) logg; c) ^; d) [Fe/H]; e) [Ca/H]; f) [Ni/H]; g) [Ba/H]; and g) 
[Eu/H]. 



the sensitivity of the 6496 A transition to departures from LTE 
appears to be a minor and secondary source of error for the M15 
target stars. 

5. Chemical Abundances in the l\115 Sample 

The wavelength range and resolution of the spectra observed for 
this study allowed for the analysis of two Ca I lines, one Ni I line, 
one Ba II line, one La II line and two Eu II Unes (see Table [2l). 
Abundances from these lines were determined using the spec- 
trum synthesis routine, synth, in MOOG. For each of the stars 
in the sample Table [TT| lists the chemical abundances and the 
associated fitting errors, which were determined as the greatest 
change in abundance for possible fits to the spectral feature. Also 
listed are the C09 parameters and chemical abundances, ^wn 
and the S/N measured here. 

In order to extract further information from the current Ml 5 
sample, groups of stars were identified with similar stellar pa- 
rameters and low Ba abundances. Due to the noise limitations, 
the key spectral features of Eu and La could not be measured 
confidently for these individual spectra. However they could be 
combined into two summed spectra such that the noise was suf- 
ficiently reduced to estimate the 'summed' heavy element abun- 
dances. Table|6]gives the mean parameters for each set of spec- 
tra including the mean [Ba/Fe] for each set, the [Ba/Fe] derived 
from the summed spectra, the upper limits on [Eu/Fe] for both 
summed spectra, and upper limit for [La/Fe] for just one of the 
summed spectra. 

Figure |5] shows the component spectra and the resulting 
summed spectra of A and B in the regions about the three key 
heavy element spectral features. The component spectra of B are 
much noisier than A as expected due to being fainter (hotter) on 
the RGB and due to it being an inherently weak line, the La spec- 
tral line profile was still too weak to be detected above the noise. 
Hence no upper limit on the La abundance could be extracted for 
summed spectrum B. The goal of the summed spectra analysis 



Table 6. Average stellar parameters of two sets of summed spec- 
tra and the resulting heavy element abundances. 





A 


B 


<S/N> 


121±25 


70±11 




4949±116 


5216+43 


<logg> 


2.1±0.3 


2.6±0.1 


([Fe/H]> 


-2. 3 ±0.1 


-2.4+0.1 


(^IV13> 


1.7 


1.6 


<[Ba/Fe]> 


-0.16±0.05 


-0.33±0.17 


[Ba/Fe] 


-0.15±0.08 


-0.39±0.15 


[Eu/Fe] 


<0.42 


<0.41 


[La/Fe] 


<0.37 




C09ID 


A= 18815, 42362, 27889 




B= 35961, 23216, 23153, 8927 



was to investigate the Eu abundance at low Ba abundances. This 
will be discussed further in the analysis of the results. 

The analysis of this sample of stars in M15 must be consid- 
ered in light of the results of the previous studies. The mean val- 
ues for the metallicity, Ba, Eu and La abundances for each key 
of Ml 5 study are given in Table |7] where the C09 metallicity 
represents the mean metallicity for this sample. 

What is noticeable is that the mean Ba abundance derived in 
this study is some 0.2 dex greater than for the majority of the 
other studies. Looking particularly at S97 and SOOb, which have 
sizeable samples likely to sample the large range in Ba abun- 
dances, there is an offset of +0.12 dex to the mean Ba abun- 
dance of S97 and an offset of +0.17 to the mean Ba abundance 
of SOOb. Both studies have dispersions in Ba abundance which 
are comparable to the spread in Ba measured for this sample. 
Hence the large range of Ba values for M15 is confirmed in this 
study but the absolute abundances appear offset to previous stud- 
ies. However it should again be noted that the 63 stars analysed 
here are the largest sample by far of M15 stars for which Ba 
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C09 IDs = 18815,42362,27889 4949±116, 2.010.2, -2.3±0.1, 121±25 




6390 6392 6493 6495 6497 6642 6644 6646 



C09 IDs = 35961,23216,23153,8927 5216±43, 2.6±0.1, -2.4±0.1, 70±11 




6390 6392 6493 6495 6497 6642 6644 6646 

WiM3lei^h(A) m\elei^h(l) WMilergth(l) 



Fig. 5. The three component spectra (blue, red, orange) and resulting summed spectrum (black) for A in the top panels, and 
the four component spectra (blue, red, orange, cyan) and resulting summed spectrum (black) for B in the bottom panels, in the 
regions of the La II, Ba II and Eu II spectral lines respectively. The C09 IDs for each set of spectra and the resulting mean stellar 
parameters {Teff, logg, [Fe/H]) and S/N for the final summed spectrum for each of A and B are given. 



Table 7. The number of stars in the sample, the mean and standard deviations for [Fe/H], [Ba/H], [La/H] and [Eu/H] for each key 
study. 







<[Fe/H]> ±0- 


<[Ba/H]> +cr 


<[La/H]> +cr 


<[Eu/H]> +cr 


W13 


63 




-2.20+0.26 


-2.05+0.15 


-1.72+0.15 


Skll* 


3 


-2.55±0.02 


-2.10+0.12 


-2.07+0.23 


-1.72+0.29 


DOlO 


57 




-2.152±0.412 






C09 


84 


-2.341±0.007 








O06 


7 


-2.40±0.05 


-2.56+0.18 


-2.46+0.25 


-1.85+0.24 


SOOb 


31 


-2.37±0.05 


-2.25+0.21 






SOOa 


3 


-2.28±0.05 


-2.28+0.19 




-1.41+0.31 


S97 


18 


-2.40±0.04 


-2.30+0.21 




-1.91+0.21 


* RGB stars only 











abundances have been determined, the next largest sample for 
Ba abundances being SOOb with 30 stars. This sample is also 
complete along the RGB based on the selection criteria for C09. 

The comparison to DO 10 is unique in that both DO 10 and 
this study are based on slightly diff'erent subsamples of C09. 
Interestingly the mean Ba abundances of these two studies are 
in reasonable agreement although the spread reported in DO 10 is 
greater As was discussed in Section|3] the ^ values determined in 
C09 resulted in a trend of Ba abundance with ^. Hence the mean 
Ba abundance from DO 10 must be considered in light of this, as 
no adjustment to the ^ values was made. The equivalent width of 
this alternate subsample of C09 were provided (D'Orazi, private 
communication) in order to make a better comparison of the Ba 
abundances using the log g - ^ relation derived here. 

Considering the mean abundances for La and Eu given in 
Table |7] the mean La abundance agrees well with Ski 1, but O06 
is overestimated in comparison. On the other hand the mean Eu 



abundance for this study lies within the range of values deter- 
mined for the other studies. 



The offset in Ba abundance between the studies can be at- 
tributable to differences in the measurement techniques and anal- 
ysis methods used to derive the stellar parameters and abun- 
dances. As a key example, the sample of 18 stars in S97 were 
analysed by the measurement of equivalents widths which were 
then used in a curve-of-growth analysis to determined the stel- 
lar parameters and abundances, whereas for C09 the stellar pa- 
rameters were determined via photometry and fiducial relations 
configured for M15. In particular, the Ba abundance determined 
in S97 was based on three Ba lines for which equivalent widths 
were measured. Two of the lines (6141.73 A, 6496.91 A) were 
typically measured to have equivalent widths at ~ 130 mA, while 
the third (5853.69 A) at -80 mA. 
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Fig. 6. C09 stellar parameters and W13 ^ attributed to the S97 
and DOlO samples, a) ^ against \ogg for C09 (black), S97ivi3 
(red) and DOlOwn (blue) and DOlO (cyan), b) HR diagram for 
C09, S97H.13 and DOlOwn- 



There is a similar difference in analysis method with DOlO 
for which equivalent width techniques were also employed in the 
analysis of the single Ba line at 6141.73 A. 

As S97 reports the equivalent widths that were measured for 
the 18 stars, and the equivalent widths from the DOlO analysis 
were also provided privately by the authors, we pursued extend- 
ing the current analysis to include the stars from these papers. 

As outlined in Table |4] six of the stars in S97 are already 
part of the current sample. Four of these stars are also part of 
the DOlO sample. In total there are 38 stars in common be- 
tween DOlO and W13. Using the corresponding C09 parame- 
ters with the modified ^ value, the Ba (6141.73 A, 6496.91 A) 
and Eu (6645 A) abundances were re-determined using the S97 
equivalent widths and the DOlO equivalent widths for the Ba 
(6141.73 A) abundance. 

Simple low order polynomial relations were derived for T^f^ 
and log g as functions of V based on the entire C09 sample. 
These functions were then used to calculate the corresponding 
Tgff and log g, and from log g the ^, for each of the 12 S97 stars 
and 19 DOlO stars. The re-analysis of the DOlO and S97 sam- 
ples using C09/W13 parameters and the resulting abundances 
are referred to as 8975^13 and DOlOiyn for the remainder of the 
paper 

Figure |6^ shows ^ compared with logg for the C09 (black), 
DO10vvi3 (blue) and S97iyi3 (red). The original ^ values used 
in DOlO (cyan) are shown for comparison. Figure |6j) is corre- 
sponding HR diagram for C09, DOlO^o and S97vyi3. 

Using these parameters, and assuming a metallicity of 
[Fe/H] =-2.341 dex (C09), stellar models were interpolated us- 
ing ATLAS9 as for the stellar models for the W13 sample. Then, 
using the equivalents widths from each study, the abundances 
for Ba and Eu were derived using the blends routine in MOOG 
taking account of hyperfine structure using the same linelists as 
used in the spectrum synthesis. 

The cross-over stars between the three studies were inves- 
tigated to look for potential systematic trends between the sam- 
ples. Table[8]lists the six stars in common between W 13, S97 and 
DOlO with their respective IDs, the Ba equivalent widths from 
S97 and DOlO, and the resulting derived abundances where pos- 
sible (8971^13, DOlOvvn). While two Eu lines were used in the 
spectrum synthesis analysis in this paper only the abundance de- 
termined for the Eu line in common with S97 (6645 A) was used 
to determine the bias between the techniques. 

Two effects can be seen here and are illustrated in Figure |7] 
which shows the A[Ba/H] between the cross-matched stars of 
W13, S97vt'i3 and DOlOi^o. Although for just a sample of six 



X 0.2 



* ★ 



* 

13.0 



^ W13(6496)-S97|^,j(6496) = 0.14+0.11 

S97^,3(6496)-S97^,3(6141)=0.14±0.34 
. W1 3(6496)-D01 (61 41 )=0.32+0.20 



15.0 
V (mag) 



★ 

17.0 



Fig. 7. Comparison of the [Ba/H] derived for the stars in com- 
mon between W13, S97h'i3 and DOlO^n as per the legend 
against V magnitude. Biases and dispersions between each sam- 
ple are listed. 



stars there is a bias of ~0. 14 dex between the [Ba/H] abundances 
derived from the specturm synthesis in W13 compared with the 
S97iyi3 [Ba/H] abundances derived from the equivalent width 
measurements. This first effects shows that even for the same 
lines differences in measurements and analysis techniques lead 
to a difference in derived abundances. 

The second effect is illustrated by the comparison of the 
S97iyi3 [Ba/H] values derived from the 6496 A line compared 
with the 6141 A line. There is also a bias of ~0.14 dex com- 
bined with a significant scatter (~0.34). Although great care was 
taken here deriving the ^ values this discrepency between the two 
strong lines measured in the same study may show the derived 
relation is still not optimised. 

Hence not only differences between measurement techniques 
and analysis must be taken account of, but also differences in the 
abundances derived by the same method but of two different Ba 
lines. This is reflected in the comparison of W13 to the DOlOvrn 
[Ba/H], which is a comparison of spectrum synthesis to equiva- 
lent width determinations as well as a comparison between the 
6496 A and 6141 A lines. The bias between these two samples is 
~0. 32+0. 20 dex, which equates to the combination of the biases 
produced by these two effects. 

In addition, while there is no particular trend of A[Ba/H] with 
V, the brighter stars show a clear offset but the fainter stars have 
a much greater dispersion and less systematic offset. At the faint 
end the noise is greater, particularly shown by the greater fit- 
ting errors measured for W13, which can for example affect the 
placement of the continuum and so the derived abundance. 

Also there are the issues between measurement and analysis 
techniques, and that it is a comparison between two different 
Ba lines. The comparison between the S97 measurements shows 
also a large scatter in values. As much as great care was taken 
in deriving the ^ values this discrepency between the two strong 
lines measured in the same study may show the derived relation 
is still not optimised. 

The bias for Eu calculated between W13 and 8975^13 is 
A[Eu/H]=0.01+0.09 dex, which is negligible. The dispersion in 
the bias reflects the random errors as expected. This is not un- 
expected as the differences in the Ba abundances are most likely 
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Table 8. Ba equivalent widths (EW) measurements from S97 and DOlO and the corresponding [Ba/H] and [Eu/H] derived from 
these measurements based on the C09 parameters and W13 ^ for the cross-over stars between W13, S97 and DOIO. The W13 
abundances are based on spectrum synthesis. 







Ba EW (mA) 




[Ba/H] 








[Eu/H] 




C091D 


S97 ID 


S97 6496 


S97 6141 


DO10 6141 


W13 6496 


S97,ri3 6496 


S91wu 6141 


DO10,ri3 6141 


W13 6645 


897^13 6645 


43788 


K825 


136 


137 


137 


-2.27 


-2.52 


-2.63 


-2.63 


-1.85 


-1.84 


40825 


K386 


129 


133 


140 


-2.42 


-2.65 


-2.72 


-2.61 


-1.94 


-1.87 


4099 


K341 


146 


146 




-2.27 


-2.42 


-2.53 




-1.67 


-1.66 


31914 


K757 


151 


147 


154 


-2.04 


-2.23 


-2.28 


-2.39 


-1.67 


-1.72 


41287 


K702 


130 


131 




-2.36 


-2.46 


-2.26 




-1.90 


-2.07 


3137 


K387 


143 


139 


133 


-2.02 


-1.96 


-2.17 


-2.28 


-1.73 


-1.65 



r.o 




-2.6 -2.4 -2.2 -2.0 -1.8 -2.6 -2.4 -2.2 -2.0 -1.8 -2.6 -2.4 -2.2 -2.0 -1.8 -2.6 -2.4 -2.2 -2.0 -1.8 

[Fel/H] [Fel/H] [Fe l/H] [Fe ll/H] 

Fig. 8. Comparisons between chemical abundances measured here with Fe from C09 for this sample: a) [Ca I/Fe] vs [Fe l/H]; b)[Ni 
I/Fe] vs [Fe l/H]; c) [Ba II/H] vs [Fe l/H]; and d) [Ba ll/H] against [Fe II/H]. Errors forCa, Ni and Ba reflect the fitting errors of the 
spectrum synthesis. Mean values and uncertainties for each sample are shown in red. 



driven by the microturbulence, the effects of which are negligi- 
ble on the weak spectral lines of Eu. The sensitivity of strong 
Ba features to microturbulence makes inter-study comparisons 
more problematic, even after considerations of differences in 
measurement and analysis techniques. The spread in the val- 
ues of A[Ba/H] are further evidence of the difficulty in obtaining 
consistency between studies. 

The goal of this exploration was to see how to combine 
these samples. The offsets in themselves could be applied but the 
large scatter, particular between W13 and DOlOvi/n for the faint 
stars, may reflect too many issues with the difference between 
the samples and techniques. It was decided that combining these 
samples potentially causes more problems then it solves but the 
samples can be analysed alongside each other for comparison as 
three distinct internally consistent analyses. 

Table |9] lists the 12 stars from S97 and the 19 stars from 
DOIO, the stellar parameters calibrated to C09 that were used to 
generate the atmospheric models and the resulting [Ba/H] (from 
the 6496 A line only for S97) and [Eu/H] abundances (for S97 
only). 

5.1. Abundance Distributions 

Ca abundances were determined for 62 of the 63 stars while Ni 
abundances were determined for 40 of the 63 stars. From the 
analysis in C09 23 of the 63 stars have derived O abundances 
and 49 have derived Na abundances (see Table [T]i. Figure |8] a. 
and b. show the distribution of [Ca/Fe] and [Ni/Fe] respectively 
compared to [Fe/H] for the M15 sample along with the sample 
mean values and standard deviations in red. 

Both Ca and Ni show a very tight distribution similar in mag- 
nitude to that of [Fe/H], implying that there is little star-to-star 
variation within the cluster of the abundances of these two ele- 



ments. As described in full in C09, O and Na both have a large 
spread of abundance values between stars in M15. However, us- 
ing the values from C09, no trend was found between Na nor O 
with any of the three heavy elements. 

Figure [8};. and d. compares the derived [Ba/H] values with 
the iron abundances determined from Fe I and Fe II lines in C09. 
There is no trend of Ba with Fe for either species and the larger 
spread of Ba compared with Fe is distinct. 

5.2. Heavy Eiements 

The primary goal of this study was to investigate Ba and Eu 
abundances in a large sample of M15 giant stars. While Ba could 
be measured for each of the 63 stars, Eu was measured for only 
20 of the stars and La was measured for just 13 of the stars. 
Figure |9] compares the derived [Ba/H], [La/H] and [Eu/H] val- 
ues with the stellar parameters. 

As discussed in Section [3] the Ba II spectral feature consid- 
ered here is very sensitive to ^ but there is clearly no trend of 
[Ba/H] with ^w\3 despite the large spread in [Ba/H] values. This, 
as well as the uniform Ca abundances, for which one of the lines 
was also sensitive to ^, confirms the appropriate derivation of the 
f-log g relation carried out for this sample. 

The observations suffered from greater levels of noise in the 
lower part of the RGB (fainter stars) hence greater uncertainties 
are reported. The subsample for which La and Eu abundances 
could be determined from their respective weak spectral features 
were necessarily located high on the RGB (brightest stars). 

Figure [10] combines the samples from the key studies of 
heavy elements in M15, including the W13 sample. The S97 
sample is that as reported in that paper Despite the range in mea- 
surement types and analysis techniques there is a clear trend of 
[Ba/H] with [Eu/H], and [La/H] with [Eu/H] even with the dis- 
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Fig. 9. Comparison of the [Ba/H], [La/H] and [Eu/H] values with the stellar parameters from C09. Errors on [Ba/H] and [La/H] are 
spectrum synthesis fitting errors. Errors on [Eu/H] are the spectrum synthesis fitting errors if just one line was measured, otherwise 
the standard deviation of the two Eu lines summed in quadrature with the fitting error. 
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Fig. 10. a) [Ba/H] and b) [La/H] compared with [Eu/H] for the 
W13 (black) sample compared with S97 (blue), SOOa (green), 
O06 (cyan), and Skll (red:RGB, pink:RHB) samples where 
available. 



persion between the studies. However, as discussed in Section|5j 
there are distinct offsets between the studies which in a basic 
comparison, such as in Figure [TO] may blur any more refined 
trend in the GC itsself. Hence for the rest of this paper only the 
re-analysis of the S97 and DO 10 samples are considered and are 
displayed separately to the W13 sample. 

Figure [TTJi & b compai-es [La/H] with [Ba/H] and [La/H] 
with [Eu/H] for the subsample of 13 W13 stars for which all 
three species could be measured, and Figure [TTt compares 
[Ba/H] with [Eu/H] for the 20 stars with both Eu and Ba abun- 
dances. Figure fTTB compares [Ba/H] with [Eu/H] for the 18 stars 
in S97. Included on each graph is the line of best fit to each 
dataset where the weighting of each point was defined as the in- 
verse of the errors summed in quadrature. The slope (m), correla- 
tion coefficient (r) and p-value (p) for each comparison is given. 



Also included in red are upper limits in the cases for which a line 
could not be accurately measured but was still detectable. The 
limits derived from the two summed spectra are also included in 
grey. 

For each comparison of the W13 samples the calculated p- 
value is less than 0.05 which indicates that the linear relation of 
these data is unlikely to have been produced by chance and is 
therefore significant in each case (i.e. the null hypothesis, that 
there is no correlation, is rejected). However the S97iyi3 sample 
fails this test. The strongest correlation is between La and Ba, 
with the La to Eu correlation being of a similar magnitude and 
both with values of the correlation co-efficient (r) being greater 
than 0.7, which indicates a strong positive linear relationship. 

For W13 the correlation between Ba and Eu is also signifi- 
cant based on the p-value, but the value of r indicates only a mod- 
erate positive linear relationship. The graph itsself shows a great 
deal of dispersion with at least three extrema, unlike the compa- 
rably clean (albeit smaller sample) relation of La to Eu. There is 
some dispersion in the La and Ba relation. The correlations of La 
to Eu have been confirmed for three other GCs, as well as Ml 5, 
in the investigation of literature sources in lRoedererl(l201 Ih . 

However for both samples, there does appear to be evidence 
of a bimodal distribution in [Ba/H] as seen in Figure [TTt & d. For 
W13, at [Ba/H] — 2.15 there is a clear separation in the distribu- 
tion, with the low Ba abundance group being more dispersed 
in both Ba and Eu than the high Ba abundance group, which 
appears to have a much tighter distribution (ignoring the high 
Ba extrema). For S97ivi3 the separation occurs at [Ba/H]~-2.30 
with the low Ba abundance group showing a more scatter though 
uniform distribution while the high Ba abundance group has a 
tight core and is more dispersed in extrema. But the W13 and 
S97iyi3 samples are in reasonable agreement, if globally offset. 
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Fig. 11. Comparisons of [La/H] with [Ba/H] and [La/H] with [Eu/H] for the 13 W13 stars with La, Eu and Ba measurements (a & 
b), and [Ba/H] with [Eu/H] for the 20 W13 stars with Ba and Eu measurements (c). The errorbars for La and Ba are the fitting errors, 
while for Eu the errorbars are the standard deviation if two lines were measured, otherwise the errorbars are the fitting errors. The 
line of best fit (grey dash) and its slope (m), as well as the correlation coefficient (r) and p-value (p) of each dataset is given. Red 
arrows are upper limits on the abundances for particular stars. Grey arrows are limits determined for the two sets of summed spectra, 
d) [Ba/H] vs [Eu/H] for the 18 S97vi.i3 stars. 



Table 9. [Ba/H] and [Eu/H] derived for the 12 stars of S97 and 
the 19 stars of DOlO that are not part of the W13 sample. For 
S97iyi3 only ^(Ba) = 6496 A and /l(Eu) = 6645 A are used. The 
C09 stellar atmosphere parameters (with [M/H] =-2.341 dex for 
897(^13) and ^v^o are listed. 



S97ID 


V 


^eff 




[M/H] 




[Ba/H] 


[Eu/H] 


K479 


12.68 


4288 


0.54 


-2.341 


2.2 


-1.99 


-1.84 


K634 


12.81 


4321 


0.63 


-2.341 


2.2 


-2.56 


-1.96 


K583 


12.83 


4327 


0.64 


-2.341 


2.2 


-2.55 


-2.12 


K490 


12.84 


4329 


0.65 


-2.341 


2.2 


-2.87 


-2.06 


K853 


12.88 


4339 


0.67 


-2.341 


2.2 


-2.65 


-1.99 


K462 


12.90 


4344 


0.69 


-2.341 


2.2 


-2.31 


-1.49 


K431 


13.03 


4377 


0.77 


-2.341 


2.1 


-2.46 


-1.95 


K144 


13.06 


4385 


0.79 


-2.341 


2.1 


-2.26 


-1.70 


K1040 


13.40 


4471 


1.01 


-2.341 


2.1 


-2.19 


-1.68 


K969 


13.45 


4484 


1.04 


-2.341 


2.1 


-1.53 


-1.95 


K169 


13.47 


4489 


1.06 


-2.341 


2.1 


-2.26 


-1.69 


K146 


13.57 


4514 


1.12 


-2.341 


2.0 


-2.53 


-2.07 


C09ID 
















34332 


13.32 


4451 


0.96 


-2.41 


2.1 


-2.42 




14574 


13.58 


4518 


1.09 


-2.34 


2.1 


-2.65 




25247 


14.00 


4623 


1.34 


-2.34 


2.0 


-2.35 




34561 


14.07 


4641 


1.39 


-2.27 


2.0 


-2.15 




29401 


14.09 


4646 


141 


-2.34 


2.0 


-2.78 




37415 


14.19 


4671 


1.43 


-2.28 


1.9 


-2.17 




9753 


14.29 


4697 


1.52 


-2.40 


1.9 


-2.73 




18913 


14.44 


4735 


1.61 


-2.46 


1.9 


-3.06 




31313 


14.72 


4805 


1.75 


-2.28 


1.8 


-2.08 




33507 


15.04 


4887 


1.94 


-2.33 


1.8 


-2.21 




41437 


15.37 


4971 


2.07 


-2.29 


1.8 


-2.68 




31956 


15.54 


5015 


2.16 


-2.34 


1.7 


-2.35 




30535 


15.58 


5024 


2.18 


-2.24 


1.7 


-2.07 




41279 


15.82 


5087 


2.35 


-2.33 


1.7 


-2.15 




28350 


15.90 


5105 


2.37 


-2.37 


1.7 


-2.86 




38382 


16.19 


5180 


2.53 


-2.20 


1.6 


-2.20 




24706 


16.40 


5232 


2.67 


-2.42 


1.6 


-3.08 




36971 


16.75 


5321 


2.76 


-2.26 


1.5 


-2.39 




20994 


16.82 


5338 


2.79 


-2.28 


1.5 


-2.51 





with the bimodal distribution. The upper limits on the two sets of 
summed spectra are also shown in grey. The two sets of spectra 
were selected as having components of low Ba abundances in 
order to explore the low Eu regime, the key line of which was 
not resolved in the majority of the individual W13 spectra. One 
set (with higher S/N) clearly falls within the low Ba abundance 
mode. The second set, comprised of much noisier spectra, lies 
well below this mode and given how noisy the spectra are should 
probably be discounted. 

5.3. Bimodal Ba distribution in M15? 

The existence of a bimodal distribution in Ba and Eu in M15 was 
raised in Sneden et al. ( 1997). A re-analysis of key sets of stars 
from ISneden et al.l 019971) iSneden et all (l2000al) . 



^ . Preston et al 

(l2006l) and lOtsuki et air (l2006l) in the study of Isobeck elal 



in the two groupings so as to emphasise rather than diminish the 
separation in [Ba/H]. 

In FigurefTTt the upper limits on the determination of Eu for 
6 noisy W13 spectra are shown in red and also fall so as to agree 



(2011) determined that such a bimodal distribution is not evi- 
dent but rather that there was a continuum of increasing Ba with 
increasing Eu. 

To investigate this further Figure [12] shows a series of 
histograms of [Ba/H] for the W13 and 00101^13 samples. 
Figure [T2h is the histogram of the W13 sample with a binsize of 
0.1 dex (blue). This is the largest M15 sample of giants star to- 
date that have been analysed for their Ba abundances, which also 
well samples along the RGB. A bimodal distribution is clearly 
seen. The histogram of the subset of bright stars from W13 (V 
< 15.2) is also shown and the bimodal distribution is still clear. 
These stars are singled out as having the lowest associated error 
and hence most accurate abundnace determination. 

Figure [T2b replicates Figure [12^ but for the DO10iyi3 sam- 
ple. This similarly large sample of 57 giant stars also well sam- 
ples along the giant branch of M15. There is no clear bimodal 
distribution but also no clear single gaussian distribution either 
The distribution appears rectangular in that all the bins, includ- 
ing the extrema, are similarly populated. That the extrema bins 
are thus populated argues against this distribution being a simple 
gaussian. Even the bright subsample for DOlOivn does not clar- 
ify the distribution. It is noted that the range of values is offset to 
W13 to more Ba-poor values as expected. 

Figure [T2b combines W13 with the further 12 stars from 
S97ivi3 and the 19 stars DOlOwn by applying the offsets spec- 
ified above which corrected the S97w\3 stars and the DOlOwn 
stars to the W13 system. The overall distribution is bimodal but 
the extrema bins are also well populated. Hence the character- 
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istics of the W13 (bimodal) and DOlO (well-populated extrema 
bins) are clear and do not erase the effects of each other However 
as the crossover samples have a noticable scatter between the Ba 
results, which implies that per star the Ba abundances from each 
line and study are not in good agreement, then combining the 
three does not produce a sample that can be considered to be 
internally consistent (See Section|5]l. 

The threshold of the bimodal distribution in the W13 lies at 
[Ba/H] — 2.20 dex separating the sample into two modes. Mode 
I being the low Ba abundance mode and Mode II being the high 
Ba abundance mode. Table [TO] lists the mean and standard de- 
viation for the three samples of W13, DOlOwn and S97tyi3. 
Statistics are also given for the two modes for W13, DOlOivn 
(assuming a threshold for DOlOiyn of [Ba/H]~-2.47 dex) and 
S97n/i3 (assuming a threshold of [Ba/H] — 2.36 dex). 

The [Fe/H] values stay relatively constant across all means 
and all samples with very low associated spread reflecting the 
tight spread obtained in C09. While the absolute mean values 
for Ba are offset from each other between the three studies all 
show the expected large spread in values. When considering the 
modes the spread of the [Ba/H] values are on a similar magnitude 
for each of the samples. 

Regarding Eu, W13 and S97h'i3 show good agreement in the 
mean for each whole sample and for the mean of Mode I. Mode 
II, the high Ba mode, shows enhanced Eu for S97vi/i3 compared 
to the mean of Mode II for W13. The spread in Eu values within 
each mode is of similar magnitude and not significantly smaller 
than that of the spread in the Ba values. While La abundances are 
only available for W13 the spread on the whole sample and for 
each mode are again similar in magnitude to the corresponding 
spread in Eu. 

Certainly for each sample and each element, there is a clear 
separation in mean values between the modes. Even for Eu and 
La the separation between the means of Mode I and Mode II are 
outside the 1 cr limit. 

5.4. s-process residuals: [Ba/Eu], [La/Eu] 

The combination of the various studies confirmed the correla- 
tion of [Ba/H] with [Eu/H], and [La/H] with [Eu/H] and the 
analysis here has confirmed a bimodal distribution in Ba and Eu 
and possibly also in La. To further explore these heavy element 
abundances in M15 Figure [13] compares [Ba/Eu] and [La/Eu] to 
[Eu/H], [Ba/H] and [La/H] for both the W13 (black) and S91wn 
(blue) samples. 

The ratios of [Ba/Eu] and [La/Eu] can be considered as the 
contribution of the i-process to the heavy element abundances 
once the contribution from the r-process (Eu) has been removed. 
For W13 in Figure [T3b and f there is no trend of either [Ba/Eu] 
nor [La/Eu] with [Eu/H]. This is not unexpected, as if all the 
contribution from the r-process has been removed then the s- 
process contribution would not trend with increasing r-process 
abundances. This is borne out by the distribution lying closer to 
the estimated sol ar r-process only yield abundance (indicated as 
red dashed lines (ISimmerer et al.ll2004l) ') than to the solar abun- 
dance (red dotted lines), which is a combination of both r- and 
s-process contributions. However the distribution does not lie ex- 
actly on the estimated solar r-process only yield implying that 
there is some contribution of the heavy elements from some 
other process. This is replicated in Figure [T3t for the S97vyi3 
sample. 

Figure [T3t is particularly interesting as it shows whether 
[Ba/Eu] shows any trend with increasing [Ba/H], thus s-process 
contribution increasing with s-process abundance. The bimodal 



distribution in Ba can be seen as a plateau for Mode I, imply- 
ing that this mode has had no further heavy element contribu- 
tion from the s-process, but a distinct linear trend of increasing 
[Ba/Eu] with [Ba/H] for Mode II. Both of these distributions for 
each mode are replicated in the S97iyi3, with the linear trend in 
Mode II appearing to extend to the outliers of each sample. This 
would imply that the stars in Mode I have undergone further Ba 
enhancement processes (i.e. the s-process). 

Is this linear relation in Mode II real? As La and Ba are 
both s-process elements then an interchange of these parameters 
should show the same linear relation. However, as La could only 
be measured for 13 of the W13 sample and was not available 
as an equivalent width measurement in S97, then the sample for 
comparison is significantly reduced. In order to aid the investiga- 
tion the 8 stars with all three Ba, La and Eu abundances that lie 
within Mode II are identified in all the W13 graphs of Figure [13] 
as magenta circles. For Figure [Till there is no convincing linear 
trend, despite the one outlier. For Figure [T3t and h, if one ex- 
cludes the high La outlier, then there is no strong argument for a 
linear relation. 

An examination of the parameters of the stars within Mode II 
showed no obvious trend of Ba abundance with ^. However for 
five of these objects there is a distinct linear trend with and 
logg. Given the almost linear relations between V and these two 
parameters, and the linear relation developed here between ^ and 
logg it is likely that these are the cause of the observed linear 
trend between [Ba/Eu] with [Ba/H]. The error analysis shows the 
sensitivity of the analysed Ba line to these parameters, which is 
accentuated for stronger/more enhanced features. Hence this is 
the most likely reason why this relation is observed in the high 
Ba Mode II, not the low Ba Mode I. 

The La line is weaker and hence not so sensitive to ^. While 
there are fewer points, the lack of a linear trend should be taken 
as confirmation that the linear trend in Ba is not real. However, 
the location of the [La/Eu] distribution with respect to the r- 
process only yields and the solar abundance does not exclude 
some contribution from the s-process or some other process that 
produces heavy elements. 

6. Heavy and light element contributions within MIS 

The existence of two modes in Ba and Eu abundances in M15 
was first proposed in S97. For both modes the mean [Ba/Eu] 
determined in S97 were found to be almost exactly the same 
(<[Ba/Eu]> — 0.40+0.05 dex) and indicative of a pollution his- 
tory dominated by contributions from the r-process. In this study 
there is a separation in the mean [Ba/Eu] values of ~0.3 dex out- 
side the uncertainties (see Table [TOl l. with the value for Mode I 
being less than that for Mode II. In the S97vi'i3 sample similar 
values to W13 were found and so a similar difference outside 
the uncertainties. But both indicate no extra contribution from 
the s-process, therefore this study agrees with S97 that both 
modes have an enrichment history dominated by the r-process 
but that Mode I has been less enriched in the heavy elements 
than Mode II. 

However when considering the s-process residuals as repre- 
sented by [La/Eu] both modes have values in close agreement 
and well within the uncertainties, with Mode I being slightly 
less enhanced. Given the greater consistency that can be obtained 
from La measurements this can be well relied upon, but still sup- 
ports the above conclusion that both Modes are r-process domi- 
nated with a possible greater enrichment in Mode II. 

SOOb confirmed the conclusions of S97 on a larger sample 
of M15 stars for which Ba only was determined of the heavy el- 
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Fig. 12. Histograms exploring the three key samples of M15 giants in [Ba/H]: a) All 63 W13 stars (blue) and 24 bright W13 stars 
(green); b) All 57 DOlOvvn stars (blue) and 21 bright DOlOvvn stars (green); and c) the W13 sample combined with the further 12 
S97tyi3 stars and the further 19 DOlOivn stars with offsets applied. 



Table 10. Mean and standard deviations of [Fe/H], [Ba/H], [Eu/H] and [La/H] derived (where possible) for the three samples of 
W13, DOlOivn and S97vyi3 for each whole sample and separated into the respective low (I) and high (II) Ba modes. The number of 
stars for each calculation is also given. 





N 




[Fe/H] 




[Ba/H] 




[Eu/H] 




[La/H] 


[Ba/Eu] 


[La/Eu] 


W13 


63 




-2.34+0.06 


63 


-2.20+0.26 


20 


-1.70+0.13 


13 


-2.04+0.16 


-0.50+0.29 


-0.24+0.21 


DO10vyl3 


57 




-2.33+0.06 


57 


-2.47+0.31 














S91wn 


18 




-2.33+0.03 


18 


-2.36+0.28 


18 


-1.78+0.17 






-0.58+0.33 




W13 Mode I 


[Ba/H]<- 


-2.20 


-2.36+0.06 


30 


-241+0.16 


7 


-1.80+0.08 


5 


-2.19+0.13 


-0.61+0.18 


-0.39+0.15 


W13 Mode II 


[Ba/H]>- 


-2.20 


-2.33+0.05 


33 


-2.00+0.16 


13 


-1.65+0.13 


8 


-1.95+0.11 


-0.35+0.21 


-0.30+0.17 


DOIOm/13 Model 


[Ba/H]<- 


-2.47 


-2.33+0.05 


29 


-2.74+0.18 














DO10».i3 Mode II 


[Ba/H]>- 


-2.47 


-2.34+0.07 


28 


-2.22+0.16 














S97h/i3 Mode I 


[Ba/H]<- 


-2.36 


-2.34+0.01 


10 


-2.56+0.12 


10 


-1.92+0.15 






-0.64+0.19 




S97h/i3 Mode II 


[Ba/H]>- 


-2.36 


-2.33+0.05 


8 


-2.11+0.19 


8 


-1.76+0.18 






-0.35+0.26 
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[Ca/H] (dex) [Ca/H] <dex) 

Fig. 14. Comparison of light with heavy element abundances for 
the W13 sample: a) [Ba/H] against [Ca/H]; b) [Eu/H] against 
[Ca/H]. 

ements. SOOa investigated in detail several heavy element abun- 
dances for three of the S97 giants, two that lie in the proposed 
Mode I and one in Mode II. The signatures that they found were 
confirmation of r-process dominated abundances. 

The heavy elements can be paired with light element coun- 
terparts based on their common sites of nucleosynthesis in or- 
der to constrain the exact nature of the source of these chemical 
abundance patterns. In particular, the r-process and a-capture 
processes (producing Ca) can both occur in Supernova Type II. 



Similarly proton-capture processes (producing Na) occur during 
the evolution up the AGB, as does the .s-process in the thermally 
pulsing stage. 

S97 determined that neither the Ba nor Eu abundances cor- 
related with the a-element abundances. In this study no conge- 
lation was found between Ba and Ca, nor between Eu and Ca 
supporting the result found in S97 as shown in Figure [T4k and 
b. The conclusion drawn in S97, that is supported here, was that 
the mass of the star resulting in the Type II Supernova that pro- 
duced the r-process element was insufficient to simultaneously 
produce the o'-elements. The Ba and Eu abundances derived here 
fa ll within the range of values of collated data and model yields 
in IChen et al.l (l2006h . Based on this the two modes of M15 are 
both indicative of pollution by Type II supernova. As concluded 
in|Chen et al. (2006) and references therein. Supernova Type II 
associated with stars of progenitor mass between 20 M© and 
40 Mq are the main source of r-process elements in the early 
galaxy. 

With regards to the s-process and proton-capture elements, 
Na was measured for S97, SOOb, O06 and Ski 1 but no corre- 
lation between Na and either Ba or La was found. Similarly no 
correlation was found between the i-process elements measured 
here and the corresponding Na abundance measured in C09. This 
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Fig. 13. a) [Ba/Eu] against [Eu/H] for 8971^13 (blu e). Dotted red line in dicates the solar abundance. Dashed red line indicates the 
expected abundance from r-process only yields (Sim merer et al.ll2004 ). b) [Ba/Eu] against [Eu/H] for the W13 sample (black). 
Grey arrows are the upper limits derived from the set of summed spectra. Magenta circles identify those points with [Ba/H] 
greater than -2.20 dex and that have measurements for Ba, La and Eu. c) As for b) but for [Ba/Eu] against [Ba/H]. d) as for b) 
but for [Ba/Eu] against [La/H]. e) As for a) but for [Ba/Eu] against [Ba/H]. f) As forb) but for [La/Eu] against [Eu/H]. g) as for 
b) but for [La/Eu] against [Ba/H]. h) as for b) but for [La/Eu] against [La/H]. 



is seen in Figure [TSh and c where both Na and Ba are highly dis- 
persed for both the W13 and DOlOivn samples. 

Because the Na and Ba do not correlate then these elements 
were produced by different sources which polluted the medium 
from which the observed stars formed, as neither Na nor Ba can 
be produced in these observed stars themselves. Hence this lim- 
its the potential AGB polluter to those with sufficient mass to 
produce Na but not Ba (in the thermally pulsing st age) and is ev- 
idence for at least two generations of stars in Ml 5 (ICarretta et al.l 
This is borne out by comparison to theoret ical yields of 
s-proc ess element abundances from AGB models. iLugaro et al.l 
(I2OI2I) present s-process yields calculated from AGB models at 
[Fe/H]=-2.3 dex, comparable to the metallicity of M15. The ex- 
pected [Ba/Eu] abundances for a range of stellar masses all ex- 
ceed [Ba/Eu] =0.0 dex (see tables 3 and 4 of that paper). This 
is therefore in support of the conclusion here that any pollution 
from AGB in M15 contributed only to the light element abun- 
dances, not to the s-process element abundances. 

T he Na-O antic orrelation is a well-studied characteristic of 
GCs (ICarretta et al.i.2010, and references therein). It is proposed 
that the first generation of stars inherited the signature of the 
supernova nucleosynthesis whereby the O abundances are en- 
hanced to above solar values while the Na is depleted below 
solar values. During the lifetimes of the massive stars of this 
population the material is processed in hot H-buming such that 
the O is depleted and the Na is enhanced. This level of depletion 
and enhancement depends on the characteristics of the individ- 
ual stars. Hence a continuum is seen for the Na-O anticorrelation 
as different levels of contribution are made to the second gen- 
eration star-forming material depending on the contributing first 



generation star. Therefore the second generation of stars bear the 
signature of the nucleosynthetic processing that occurred within 
the first generation stars. 

Figure [TSb and e show the Na-O distribution for the W13 
and DO10ivi3 samples respectively. The Na-O anticorrelation is 
clear in both samples. Hence evidence of at least two generations 
of stars in M15. So where does the Ba bimodal distribution fit 
in? Does one of the modes represent the first generation and the 
other the second generation? If so then one mode would show the 
signature of super-solar O and sub-solar Na abundances, while 
the other is depleted in O and enhanced in Na. 

Figure [TSt and d decompose the Na-O distribution into the 
two Ba modes for W13 and DOlOvi/n respectively, as defined 
in FigurefTSh and c. In both instances we restricted the modes to 
those stars most likely to be in each mode, and so did not include 
those stars near the mean of each sample as illustrated by the red 
dashed lines in the figures. What is most clearly seen in the W13 
sample is that both Ba modes show the anti-correlation in Na-O. 
In the DO10ivi3 sample the anticorrelation is clear for Mode 1 
but more truncated in Mode 11. The two modes are less evenly 
populated for DOlOivn, as was seen in Figure [12] and so the 
sample size for Mode II may be insufficient. 

Based on the W13 sample, the polluters of Mode II appear to 
have undergone more processing of the material because the cor- 
relation extends further into the enhanced Na,depleted O regime. 

Taking the two Ba modes and then the two stellar generations 
within each mode there is a potential scenario for the formation 
of M15 which involves four populations of stars: 

- Mode I: Primordial showing SN enrichment (relatively poor 
Ba & Eu, sub-solar Na, super-solar O) 



16 



C. C. Worley et al.: Ba & Eu Abundances in M15 giant stars 



a) W12 



b) W12 



c) W12 



-1.6 
-1.8 
-2.0 
? -2.2 

IS 

ffi 

" -2.4 
-2.6 
-2.8 
-3.0 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 
[Na/Fe] 



0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 



-0.2 0.0 0.2 0.4 0.6 0.8 
[0/Fe] 



0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 



• Mode 1: 13 stars 

• Mode II: 11 stars 



• * 



-0.2 0.0 0.2 0.4 0.6 0.8 
[0/Fe] 



d) DO10 




e) DO10 



f) DO10 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 
[Na/H] 




-0.2 0.0 0.2 0.4 0.6 0.8 
[0/Fe] 



0.8 
0.6 



if 0.4 



0.2 
0.0 
-0.2 



• Mode 1: 15 stars 

• Mode II: 7 stars 



-0.2 0.0 0.2 0.4 0.6 0.8 
[0/Fe] 



Fig. 15. a) [Ba/Eu] against [Na/Fe] for 49 stars in the W13 sample (black). Dashed grey line indicates the W13 mode threshold. 
Dash-dotted grey lines indicates the mean cr[Ba/H] above and below the threshold. Blue points identify those clearly in Mode I, 
red points are those clearly in Mode II. b) [Na/Fe] against [O/Fe] for the W13 stars with both measurements, c) As for b) but the 
red and blue points from a), d) to f) are as for a) to c) but for DOlOiyn (magenta). 



- Mode I: Second generation showing SN enrichment and H- 
burning products (relatively poor Ba & Eu, enhanced Na, 
depleted O) 

- Mode II: Primordial showing SN enrichment (relatively rich 
Ba & Eu, sub-solar Na, super-solar O) 

- Mode II: Second generation showing SN enrichment and H- 
burning products (relatively rich Ba & Eu, enhanced Na, de- 
pleted O) 

In such a scenario each mode was initially polluted by the ex- 
plosive nucleosynthesis of massive stars with slightly different 
heavy element yields resulting in the bimodal distribution that 
has been observed in the W13 sample. This difference in com- 
position may have had the follow on effect of greater processing 
of the O to Na for the higher Ba mode (II). Integral to formulat- 
ing such a scenario are the timescales that have been involved. 
Did these two modes evolve simultaneously, and somehow dis- 
connectedly either separately or within the same cluster space, 
or did one population evolve followed by the other? 

Alternatively these signatures could reflect localised events. 
The stars all formed on the same timescale but the initial primor- 
dial material was not well-mixed such that stars in one region 
were more enhanced (depleted) in the heavy elements depending 
on the level of polllution by the SN. If the bimodal distribution 
is real then it would appear that the SN events were restricted 
to two types. If the bimodal distribution is not real, and there is 



a continuum of heavy element abundances, as might be argued 
based on the DOlOivn sample, then this restriction is relaxed 
and the Na-O anticorrelation may be representative of a single 
mode scenario involving a single primordial population and a 
single second generation population. 

Ski 1 investigated the heavy element abundances for 7 RHB 
and 3 RGB stars in M15, in combination with a re-analysis of 
the S97 stars using the improved analysis techniques, and found 
no evidence of a bimodal distribution in Ba and Eu for the sam- 
ple. Detailed investigation was given into light as well as heavy 
5-process elements whereby the signatures of the 10 key stars 
showed evidence of enrichment in the heavy elements that was 
not attributable to either the r- nor i-process but potentially due 
to a Lighter Element Primary Process (LEPP). Due to the wave- 
length range of the observations analysed in the current paper 
there were no spectral features of any light i-process elements 
(Sr, Y, Zr) present, hence this line of research could not be pur- 
sued. 

The above investigation of the s-process residual with 
[Ba/H], [La/H] and [Eu/H] confirmed that both modes are dom- 
inated by contributions from the r-process, hence pollution by 
explosive nucleosynthesis in massive stars. Some small contri- 
bution to the heavy element abundances by another process can- 
not be discounted as the residuals for both modes lie above the 
r-process only yield limit, although predictions of heavy element 
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yields from explosive nucleosynthesis have their own set of un- 
certainties. As such the distribution of heavy elements in this 
sample of Ml 5 giant stars implies one mode that is metal poor 
with low Ba, La and Eu abundances, and a second mode which 
is slightly less metal poor (A[Fe/H]=+0.03 dex) with relatively 
higher Ba, La and Eu abundances that all have similar disper- 
sions on the order of the dispersion for Fe. 

In effect the analysis carried out here has drawn the same 
conclusions as S97, SOOa and SOOb that there is a bimodal dis- 
tribution in the Ba, Eu, and now potentially La, abundances in 
M15 giant stars and that both modes therefore have a heavy ele- 
ment enrichment history that is dominated by the r-process. The 
addition of the Na-O abundances from C09 have added another 
layer of information which potentially describes a very complex 
formation scenario for this very interesting GC. 

Throughout Figures [TT] and [13] there is an outlier (C09 ID 
42262) in the high Ba Mode that is enhanced in Ba and La to a 
level which is significantly greater than the rest of the W13 stars. 
However, the Eu abundance is fairly typical of that mode. The 
s-process residuals for this star place it close to the solar value. 
Given the conclusions above and that because this is an RGB star 
it cannot be producing the heavy elements internally, the most 
likely scenario is that this star has been the recipient of mass 
transfer from a companion star that has resulted in the observed 
enhancments. The fact that the star is not further enhanced in Eu 
(the r-process element) and both Ba and La are heavy s-process 
elements implies that the polluting star had evolved through the 
thermally pulsing AGB stage. 

One other stai- (C09 ID 37931) in the W13 sample is simi- 
larly enhanced in Ba but no La abundance could be measured. 
This star is also enhanced in Eu, albeit with large uncertainties. 
There are potentially two such stars within the S97 sample as 
well. 



7. Conclusions 

This is the largest study of barium in giant stars in M15 to- 
date. The sample, based on the original programme in C09, is 
well-sampled along the RGB. La and Eu abundances were de- 
termined for a subset of this sample furthering this investigation 
into heavy element abundances in M15. In parallel the equiva- 
lent width measurements from DO 10 and S97 were re-analysed 
on the C09 scale which provided an interesting comparison for 
key element distributions. 

The stellar parameters for these objects were take from C09, 
except for ^ which was re-determined based on a ^-log g relation 
derived from a subsample of stars for which a reasonable sample 
of Fe I lines could be measured. The key Ba line under analysis 
was shown to be highly sensitive to ^ and so a great effort was 
made to reduce the associated errors. 

A bimodal distribution of Ba, Eu, and potentially La, was 
found in the W13 sample in agreement with previous stud- 
ies. Based on the W13 results, an examination of the r- and s- 
process contributions to these heavy element abundances leads 
to the two modes having similar pollution scenarios of a single 
burst of explosive nucleosynthese by massive stars (for exam- 
ple. Supernova Type 11) that enriched the interstellar medium 
from which the stars in each mode formed. The enrichment was 
to different degrees between the modes. 

Further enrichment by another process that produces heavy 
elements can not be discounted but the wavelength region ob- 
served here did not allow for analysis of other heavy element 
features. 



For both modes a Na-O anticorrelation was observed which 
is indicative of two stellar generations in each mode: a primodial 
and second generation in Mode 1, and a primordial and second 
generation in Mode 11. A scenario that can account for such a 
set of stellar populations is necessarily complex. If the Ba dis- 
tribution is more indicative of a continuum of Ba and Eu abun- 
dances then a single mode of stellar evolution could explain the 
observed Na-O anticorrelation. 

While this study strongly supports the argument for a bi- 
modal distribution in Ba, La and Eu for M15, there are inher- 
ent uncertainties in the measurement of the key Ba line due to 
its sensitivity to ^. La would be a better candidate as it has no 
such sensitivites. Hence re-observing the entire sample consid- 
ered here at higher resolution and S/N in order to complete the 
La and Eu abundances for the sample would prove definitive. 
Extending the observed wavelength regions to include features 
of the light i-process peak (Y, Sr, Zr) would enrich the analysis 
as to the potential contribution from other heavy element pro- 
cesses. 
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Table 11. Chemical abundances, uncertainties and number of lines measured here for the sample of 63 M15 giants stars. TheC09 stellar parameters, ^wi3 and measured 

S/N are also listed. 
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0.02 


0.10 




0.02 


0.10 
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84 
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2.47 
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0.10 


0.10 




0.01 
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ID S/N Tgff logg [Fe/Hl fvvB [Ca/Fe] cr N [Ni/Fe] cr [Ba/Fe] cr [La/Fe] cr [Eu/Fe] cr N 
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64 
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2.57 


-2.35 


1.6 


0.00 


0.13 1 


0.11 


0.25 


0.11 


0.10 
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55 
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2.65 
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- 


- 


- 


- 
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0.10 
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2.59 
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0.13 1 


- 


- 


-0.09 
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0.10 
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0.06 


0.13 1 






0.18 


0.10 


3825 


50 


5323 


2.76 


-2.26 


1.5 


0.11 


0.13 1 






0.52 


0.10 


40762 


58 


5325 


2.78 


-2.34 


1.5 


0.19 


0.13 1 






0.32 


0.13 


8329 


53 


5344 


2.86 


-2.41 


1.5 
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